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1 Introduction 

In the standard model (SM), the electroweak symmetry-breaking takes place through the Higgs 
mechanism in which a complex scalar doublet with a non-zero vacuum expectation value is in- 
troduced and the existence of one scalar particle, the Higgs boson, is predicted QUI). The gauge 
bosons derive their masses from the additional degrees of freedom gained from the symmetry 
breaking, and the fermions acquire mass through the direct interaction with the Higgs field 
itself. It is possible that the mechanism that generates the fermion masses is independent of the 
Higgs boson. Such a Higgs boson is usually referred to as fermiophobic (FP) [3 . 4j. Its decay to 
W and Z bosons proceeds as in the SM, while the decay to photons proceeds via W loops, i.e. 
decays to photons via fermion loops are excluded by the model. Since decays to bb and tt are 
forbidden at tree-level, the branching fraction for a low mass FP Higgs boson (?«h ~ 120 GeV) 
to decay to two vector bosons or two photons is enhanced by an order of magnitude with re- 
spect to the SM ISHZl. Previous searches at LEP [8], the Tevatron [9|, and the LHC [10 1 rule out 
an FP Higgs boson lighter than 121 GeV at 95% confidence level (CL). 

In this letter we report on a search for an FP Higgs boson in the mass range 110-300 GeV, 
in the 77, WW, and ZZ decay modes in proton-proton collisions at a centre-of-mass energy 
of 7TeV with data collected in 2011 by the Compact Muon Solenoid (CMS) detector at the 
LHC. The production of the FP Higgs boson is by vector boson fusion (VBF) and associated 
production with a vector boson (VH). While for the cases of WW and ZZ the search relies on a 
re-interpretation of the standard model Higgs boson searches, for the 77 final states a dedicated 
analysis has been put in place. The descriptions of the analyses emphasise the sub-channels and 
techniques not previously described in recent publications of SM analyses |TTHT6| , namely: the 
lepton tag and the use of a two-dimensional fit in the 77 decay channel, and the lepton tag in 
the WW decay channel. 

2 The CMS detector 

While the CMS detector is described in detail elsewhere ffT7| , the key components for this anal- 
ysis are summarised here. 

The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal diame- 
ter, providing a field of 3.8 T. Within the field volume are a silicon pixel and strip tracker, a crys- 
tal electromagnetic calorimeter (ECAL) and a brass /scintillator hadron calorimeter (HCAL). 
Muons are measured in gas-ionisation detectors embedded in the steel return yoke. Extensive 
forward calorimetry complements the coverage provided by the barrel and endcap detectors. 

The inner tracker measures charged particles within the pseudorapidity range \rj\ < 2.5. It 
consists of 1440 silicon pixel and 15 148 silicon strip detector modules. 

The ECAL consists of two parts: the barrel which covers the pseudorapidity t] range \ t]\ < 1.48, 
and the endcaps covering the range 1.48 < \rj\ < 3.0, where t] = — In [tan (9/2)], and 6 is the 
polar angle of the anticlockwise trajectory of a particle with respect to the beam direction. When 
used for detector positions the trajectory is assumed to originate at the nominal interaction 
point, corresponding to the coordinate system origin. The ECAL consists of lead tungstate 
crystals arranged in a quasi-projective geometry. In the barrel region the front face of the crystal 
is approximately 22 x 22 mm 2 , corresponding to a granularity of Atj x Acp = 0.0174 x 0.0174. 
In the endcap the front face of the crystals is approximately 29 x 29 mm 2 . 

In the region \t]\ < 1.74, the HCAL cells have a granularity of Arj x Acp = 0.087 x 0.087, and 
for \t] I < 1.48, map onto arrays of 5 x 5 crystals in ECAL to form calorimeter towers projecting 
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radially outwards from the nominal interaction point. At larger values of \rj\, the size of the 
towers increases and the matching ECAL arrays contain fewer crystals. Within each tower, 
the energy deposits in ECAL and HCAL cells are summed to define the calorimeter tower 
energies, subsequently used to provide the energies and directions of hadronic jets. A quartz- 
fibre Cherenkov calorimeter extends the coverage to \t]\ < 5.0. 

3 Search channels 

The search is performed using three decay modes of the FP Higgs boson. In the analysis of the 
77 and the WW decay modes, characteristic signatures of FP Higgs boson production via VBF 
and VH are exploited to select events and suppress background: the two forward jets produced 
by the scattered quarks in VBF production (dijet tag) and isolated charged leptons (electrons or 
muons) from decays of the vector bosons in VH production (lepton tag). 

In the 77 decay mode, events with two isolated and high transverse momentum (pj) photons 
are selected. Seven event classes are defined: three of diphoton events which additionally re- 
quire at least a pair of jets, an isolated muon, or an isolated electron, respectively, and four com- 
prising the remaining diphoton events which are subdivided according to the photon shower 
shape and position in the detector |fTT| . 

In the WW decay channels, the events are characterised by the presence of two opposite sign, 
isolated high pj leptons from the W decay, together with large missing transverse energy (E™ 1SS ) 
due to undetected neutrinos. One sub-channel is defined by additionally requiring the pres- 
ence of two jets with the VBF topology [12|, and another by requiring the presence of a third 
isolated, high pj lepton. The selection requiring a third lepton is intended to select Higgs boson 
production in association with a W boson (WH). Final states with one or zero jets, each sepa- 
rated into events where the leptons have the same flavour and those where there is an electron 
and a muon, add a further four WW sub-channels, giving a total of six. 

In the case of the ZZ final state, production mode signatures are not exploited, and the results 
of SM Higgs searches |fl3i - [T6Tl , comprising 19 sub-channels, are simply re-interpreted from the 
different signal rates expected in an FP Higgs boson model. 

The final result is obtained from the combination of 32 mutually exclusive sub-channels from 
the three decay modes, 77, WW, ZZ, as summarised in Table [l] The luminosity calculation 
for the datasets used has been updated with respect to that used in the SM Higgs production 
analysis of the same channels, published in the references given in the table. 

The cross sections for the Higgs boson production mechanisms and decay branching fractions, 
together with their uncertainties, are taken from Ref. [18] and are derived from Refs. |19TI24| . 
The VBF production of the FP Higgs boson signal is simulated using the next-to-leading order 
matrix element generator POWHEG 1.0 interfaced with PYTHIA 6.4.24 [251 for parton showering 
and fragmentation. The VH signal production channel is simulated with PYTHIA. Samples of 
Monte Carlo (MC) simulated events used in the analysis are passed through the Geant4 [26 1 
model of the CMS detector and reconstructed with the same software as used for collision data. 

3.1 Diphoton 77 decay mode 

In the 77 channel [TrT| , two isolated photon candidates are required to be within the ECAL 
fiducial region \t]\ < 2.5, excluding the barrel-endcap transition region 1.44 < \t]\ < 1.57. 
The shape of electromagnetic shower is used to identify photons, while track veto is used to 
exclude electrons. Isolation is used to reject the background due to electromagnetic showers 
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originating in jets - mainly due to single and multiple 7i°s |TT). The isolation requirements are 
applied as a constant fraction of the candidate photon pj, effectively cutting harder on low pj 
photons. The R9 variable, defined as the energy sum of 3x3 crystals centred on the crystal with 
maximum energy deposit divided by the total clustered energy is used to distinguish photons 
of well measured energy. 

3.1 .1 Dijet tag event class 

Candidate diphoton events for the dijet-tagged channel have the same selection requirements 
as in the SM search |TT| . In the events from the VBF production, the pj of the Higgs boson is 
boosted giving enhanced asymmetries in the photon pair energies and hence favoring a lower 
threshold on one of the two photons. The threshold requirements for this class are > 
55 x m 77 /120, and p^(2) > 25GeV. 

For each event, hadronic jets are clustered from the reconstructed particles with the infrared 
and collinear safe anti-/c t algorithm [27], operated with a size parameter R of 0.5. The selec- 
tion variables for the jets use the two highest pj jets in the event with pseudorapidity \r/\ < 
4.7. The selection requirements are optimised to obtain the best expected limit at 95% CL on 
the VBF signal cross section with fully simulated VBF signal events and the diphoton back- 
ground estimation from data |(TT][. The pj thresholds for the two jets are 30 and 20GeV, and 
the pseudorapidity separation between them is required to be greater than 3.5. The dijet mass 
is required to be greater than 350 GeV. Two selection criteria, relating the dijet to the diphoton 
system, are applied: the difference between the average pseudorapidity of the two jets and the 
pseudorapidity of the diphoton system is required to be less than 2.5 |28], and the difference 
in azimuthal angle between the diphoton system and the dijet system is required to be greater 
than 2.6 radians. 

3.1 .2 Lepton tag event classes 

Candidate diphoton events for the lepton-tagged channel have the same selection requirements 
imposed on the photons as in the SM search |HT| except for the pj thresholds. As it is the case in 
the VBF, the pj of the Higgs boson is also boosted in the VH production. To maximize the signal 
efficiency, the photon pj thresholds are set to Pj(l) > 45 x m 77 /120, and Pj{2) > 25 GeV. 

The lepton tag requires at least one muon or electron with pj > 20 GeV, within < 2.4 for 
muons, and \tj\ < 2.5 for electrons. Electrons are identified as a primary charged particle track 



Table 1: Summary of analysis channels and sub-channels included in the combination. 



Channel 


tny[ range 
(GeV) 


Sub- 
channels 


Luminosity 

(fb- 1 ) ' 


Reference 


H— > 77 


110-150 


4 


5.1 




H— > 77 + dijet 


110-150 


1 


5.1 


111] 


H— > 77 + lepton 


110-150 


2 


5.1 




H^WW^ 2£2v 


110-300 


4 


4.9 


112] 


H-fWW^ 2£2v + dijet 


110-300 


1 


4.9 


ma 


H-fWW^ 2£2v + lepton 


110-300 


1 


4.9 






110-300 


3 


5.0 


m 


H-rZZ-? 2£2v 


250-300 


2 


5.0 


EI 


H^ZZ^ 2£2q 


130-165, 200-300 


6 


5.0 


ma 


H-^ZZ^ 2£2t 


180-300 


8 


5.0 


116] 
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and one or more ECAL energy clusters corresponding to this track extrapolation to the ECAL 
and to possible bremsstrahlung photons emitted along the way through the tracker material. 
Muons are identified as a track in the central tracker consistent with either a track or several 
hits in the muon system, not associated with a significant energy deposit in the calorimeters. 

The leptons are required to be isolated, using isolation criteria similar to those used for pho- 
tons [11], and to be separated from the photons by AR > 1, where AR = y 7 At] 2 + Acp 2 . To 
protect against the background events that arise from an electron misidentified as a photon in 
the SM process Z — > ee, it is required that the mass of the photon-electron system is not within 
±5 GeV of the nominal Z boson mass. 

3.1 .3 Untagged event classes 

A substantial fraction of the H — > 77 signal events are not expected to pass either the dijet 
or lepton tag. A statistically independent search is performed on untagged events by using 
diphoton events that pass the 77 selection, and photon pj requirements of Pj(l) > m 77 /3 and 
Pj(2) > m 77 /4, but do not pass the selection for either of the two tagged channels. Higgs 
bosons produced by VBF and VH mechanism have a harder transverse momentum spectrum 
than those of the photon pairs produced by the background processes II29I , and thus the back- 
ground can be rejected while retaining high signal efficiency, by placing a requirement on the 
transverse momentum of the diphoton pair, pj 7 . It is required that 7ij 7 = Pj 7 / m 77 > 0.1. 

The selected events are divided into four classes according to the expected mass resolution and 
amount of background contamination |TT1| . Two photon classifiers are used: the minimum Rg 
of the two photons, RJf m , and the maximum absolute pseudorapidity of the two photons. The 
class boundary values for Rg and pseudorapidity are chosen to match those used to categorise 
photon candidates for photon identification cuts. The untagged diphoton event classes are: (a) 
both photons in barrel and R9™ > 0.94, (b) both photons in barrel and R™ in < 0.94, (c) one or 
both photons in endcap and R™ m > 0.94, and (d) one or both photons in endcap and R™ in < 
0.94. 

Table 2: Number of selected events in the 77 event classes, for data in the mass range 
100 — 180 GeV and for an FP Higgs boson signal (?hh = 120 GeV). The expected number of 
background events in the signal region 115-125 GeV obtained from the background fit and the 
mass resolution for the 120 GeV FP Higgs boson signal in each event class are also given. 





Dijet 
tag 
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Unta 


gged 




(a) 


(b) 


(c) 


(d) 


Data (100 < m 77 < 180 GeV) 


122 


9 


3866 


5496 


3043 


4201 


Signal (m H = 120 GeV) 


21.8 


4.4 


23.1 


23.9 


10.1 


11.5 


Expected bkg (115 < m 77 < 125 GeV) 


21.2 


1.3 


678.5 


985.2 


537.5 


754.3 


a ef t (GeV) 


1.67 


1.63 


1.19 


1.70 


2.54 


2.94 



The numbers of events in the 77 sub-channels are shown in Table |2j for simulated signal events 
and for data. A Higgs boson with Mh = 120 GeV is chosen for the signal, and the data are 
counted in the mass range 100-180 GeV. The table also shows the mass resolution, t7 eff , defined 
as half the width of the narrowest window containing 68.3% of the distribution. 

3.1.4 Signal and background modelling 

In the SM H — > 77 analysis [11], the diphoton mass spectrum of the signal and background 
are assumed to be described by analytical functions. The signal shape was determined from 
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the Z to electron mass spectrum, while the background was described by smoothly falling 
analytical functions of various forms. In this analysis an additional observable based on the 
diphoton transverse momentum, 7ij 7 , is used to construct two-dimensional (2D) PDF for the 
four untagged event classes. This enables further exploitation of the difference between the 
signal and background diphoton transverse momentum spectra. Z boson decays to electrons 
are used to derive the amount of additional smearing that needs to be applied to photons in 
MC simulated events to reproduce the energy resolution observed in the data. These smearing 
corrections are between 0.7 and 3%, derived for photons separated into four rj regions (two 
in the barrel and two in the endcap) and two categories of R9. The uncertainties on these 
corrections and factors accounting for the difference between photons and electrons are taken 
as systematic errors in the limit setting procedure. 

The signal mass model, A4 s (m 77 ), is extracted, after the smearing, by parameterising the m 77 
distribution in simulated signal events with a sum of Crystal Ball [30J and Gaussian functions. 

In the untagged 77 event classes, where a 2D analysis using m 77 and TLp is performed, it is 
necessary to define a signal model that is a function of these two observables, in the regions 
100 < m 77 < 180 GeV and 0.1 < TLp < 2. The correlations between the two variables are 
neglected, because the mass resolution of the Higgs boson has little dependence on its momen- 
tum. The probability distribution function (PDF) becomes a product of two one-dimensional 
PDFs, one for mass, M. s {m 77 ), and one for the second observable, /C s (7Tj 7 ), empirically de- 
rived to be a sum of Gaussian and bifurcated Gaussian (a Gaussian PDF with different widths 
on left and right side of maximum value) functions. The nj shape uncertainty contributes less 
than 1% to the expected exclusion limit. 

3.1 .5 Background modelling in dijet and lepton tag classes 

For the dijet-tag event class, the background model is derived from data, by fitting the diphoton 
mass distributions over the range 100 < m 77 < 180 GeV. The choice of the function used to 
fit the background and the choice of the range are made based on a study of the possible bias 
introduced by the two choices. The bias is studied for both the limit in the case of no signal and 
the measured signal strength in the case of a signal. 

In case of the lepton-tag class, the requirement of an additional isolated lepton suppresses the 
contribution from QCD background processes. The remaining background is small, coming 
predominantly from electroweak processes. Its shape is derived from fitting the MC simulation 
for muon tags, and from fitting a combination of data and MC simulated events for electron 
tags. For the electron-tagged events, a control sample (CS) is derived from data by requiring 
one of the photons to be matched with a track. This CS represents the reducible background 
with enhanced statistics. The final shape for this channel is the sum of the fits to simulation 
and the CS, with the two components weighted by the cross sections of the main irreducible 
and reducible background processes. The sum is normalised to the data yield in the range 
100-160 GeV. 

Bias studies are performed using a number of generated pseudo-experiments with background 
only and signal plus background hypotheses. It is observed that using a second order polyno- 
mial fit to the range 100 < m 77 < 180 GeV for the dijet-tagged events results in only a small bias 
in either excluding or finding a Higgs-boson signal in the mass range 110 < m 77 < 150 GeV. 
In both cases the maximum bias is found to be at least five times smaller than the statistical 
uncertainties of the fit. For the lepton-tagged classes, it is found, from the same technique, that 
the small number of selected events allow the use of exponential functions for the fits without 
introducing any significant bias as compared to the statistical uncertainty. 
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The data distribution of m 77 and the corresponding background models in the three tagged 
event classes are shown in Fig. [TJ For the dijet-tagged class the statistical uncertainty bands 
computed from the fit are shown. For both the muon-tagged and electron-tagged classes, un- 
certainty bands are not shown. For these two classes the dominant statistical uncertainty on 
the background model is obtained from the number of events which are used to determine the 
normalisation of the fit. 
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Figure 1: Data and the background model fits to the m 77 distribution for the three diphoton 
event classes from tagging channels: for dijet-tag with a second order polynomial fit (top), 
muon-tag (middle) and electron-tag (bottom) with an exponential fit for the lepton-tag. For the 
two bottom plots the uncertainty bands are not shown, as explained in the text. Signal for an 
FP Higgs with a mass of 120 GeV is overlaid for reference. 



3.1 .6 Background modelling in the untagged classes 

In the untagged 77 event classes, where a 2D analysis using m 77 and rCf is performed, the 
background model is a distribution of these two observables. The nominal background PDF 
accounts for a linear correlation between the two observables and has the following form: 

Vb = Mb{m 71 , 7r^ 7 |fl /«i) x fcb{nj 7 ) = m 77 +<!l7rT x fCbirCj 1 ) (1) 
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The empirical background PDF for the second observable, ^(tZj 7 ) is defined as a sum of an 
exponential function {£) of slope Tg and a Gaussian function (Q) of width ctq and mean fixed at 
zero: 

1C b (nP\T Br f d ,cr G ) = f d £{n?\T B ) + {\-f d )G{n^\0,a G ) (2) 

A power law function is chosen to describe the m 77 distribution. The data in each of the four 
untagged event classes are fitted separately. Figures [2] and [3] show the data and the fit results 
projected on m 77 , and on TZj 7 , respectively, for each class. 

Goodness-of-fit tests are performed measuring the bias of the model due to correlations be- 
tween the two observables and due to choice of functional forms. Pseudo-experiments are 
performed, generated from alternative background models, and the signal plus background 
model is fitted for various test masses. The bias is taken as the mean of the pull distribution, 
which is defined as the difference between the fitted and generated signal strength divided by 
the statistical error from the fit in each event. If the pseudo-experiments are generated from a 
background model containing a linear correlation of the two variables, a maximum bias of 60% 
is observed if a linear correlation is not included in the fitting function. Tests have shown that 
with a linear plus quadratic correlation in the model a fit with only a linear correlation results 
in a bias of less than 13% in the entire fit range. This bias is regarded as negligible and thus it is 
concluded that a fit function with linear correlation is adequate. With a similar procedure, the 
nominal fit function is tested against alternative models of mass shape with functions including 
linear correlations. A maximum bias of about one quarter of the statistical error is measured, 
which is negligible. 

3.2 Diboson WW decay mode 
3.2.1 Dijet tag event class 

The H — > WW**) — > 2£2v analysis [12J selects events with two isolated leptons of opposite 
charge, large missing transverse energy, and two jets with VBF topology. One lepton is required 
to have pj > 20 GeV, while the second is required to have pi > 10 GeV. The fiducial region is 
1 1] | < 2.4 for muons and 1 1] | < 2.5 for electrons. If the two leptons have the same flavour, the one 
with lower pj is required to have a pj of at least 15 GeV to suppress the Drell-Yan background. 
The missing transverse energy requirement is applied by means of a selection on the projected 
£miss defied as the component of £™ ss transverse to the nearest lepton if that lepton is within 
n/2 in azimuthal angle, or the full E™ 1SS otherwise |31], which is required to be larger than 
about 40 GeV, the precise value depending on the number of vertices found in each event. Jets 
are required to have a pj > 30 GeV and < 5, and the two jets with the highest pj are 
chosen as tag jets ]•£). The VBF selection requires — rjj 2 \ > 3.5 and tn^h > 450 GeV. The 
distributions of and Arjjj = \rjj l — t]j 2 \ are shown in Fig. Rafter the WW selection. Besides 
the cuts on the dijet system, it is required that the event has no other jets above 30 GeV found 
between the tag jets in pseudorapidity. These selections suppress the dominant background 
coming from top-quark production, which is also reduced by requiring a b-veto on the tag jets 
using jet impact parameters. The Drell-Yan background is suppressed by requirements on the 
dilepton system mass and momentum, as well as on the angle between the dilepton and the 
dijet system. Events with additional leptons above 10 GeV are rejected. After all requirements, 
between 10 and 20 data events remain in the signal regions which are defined according to 
the Higgs boson mass hypothesis fi"2 | . These events are compatible with a background only 
hypothesis. The main backgrounds, namely tt and Z + jets, are estimated from the data: the 
tt is measured in the region where the central jet is b-tagged, while the Z + jets is measured 
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Figure 2: m 77 distribution in data (points) in the four categories of the untagged 77 sub- 
channel (a)-(d) defined in Section |3.1.3 together with background model fits of a power 
function including linear correlation to Ttp '. An MC simulated FP Higgs boson signal (oth= 
120 GeV) is overlaid for reference. 
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Figure 3: Ttp distribution in data (points) in the four categories of the untagged 77 sub- 
channel (a)-(d) defined in Section 3.1.3 together with background model fits of a sum of expo- 
nential and Gaussian functions centred at zero. A MC simulated FP Higgs boson signal (oth= 
120 GeV) is also shown, scaled by 200 and modelled with a sum of a Gaussian and a bifurcated 
Gaussian. 
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Figure 4: The distributions of m /1/2 (left) and At]jj = \ rj^ — m % | (right) after the WW selection. 
Contributions from backgrounds are shown together with data, and the expected signal from 
VBF production is shown for m\\ = 160 GeV. The arrows show the minimum values in the 
selections applied for the WW dijet-tag analysis. 



under the Z peak. Also the contamination coming from W+jets and QCD is measured in data, 
in a phase space with relaxed lepton identification. The WW background is evaluated from 
simulation. 

3.2.2 Zero and one-jet event classes 

The results obtained by the inclusive analysis |fl2j in the zero and one-jet classes are also in- 
cluded in the limit extraction. The electron and muon selection is the same as for the dijet- 
tagged class. A pj threshold of 30 GeV is placed on the jets, the number of which defines the 
classes. Events are split into same-flavour and opposite-flavour dilepton sub-channels, since 
the background from Drell-Yan production is much larger for the same-flavour dilepton events. 
The dominant background for these classes is from WW, together with Z + jets and top-quark 
production in the one-jet class, as well as contaminations from W + jets and QCD. 

3.2.3 Lepton tag event class 

The WH — > WWW — > 3l3v analysis selects events with three charged leptons, either electrons 
or muons, large £™ 1SS , and low hadronic activity. The third lepton has to be isolated and have 
pj > 10 GeV, and it is required that there be no jet with pj > 40 GeV in the event. The dominant 
background comes from WZ — > 3£v production, which is largely eliminated by requiring that 
the invariant mass of all same-flavour oppositely charged lepton pairs is not within ±25 GeV of 
the nominal Z boson mass. In addition, the smallest dilepton mass constructed from oppo- 
sitely charged leptons is required to lie between 12 and 100 GeV, and the smallest distance, AR, 
between them is required to be less than 2. The background processes with jets misidentified as 
leptons, e.g. Z+jets and top, as well as the WZ — > 3£v background are estimated from data. The 
small contribution from the ZZ — > 41 process with one unreconstructed lepton is estimated us- 
ing simulated samples. After all cuts, 7 data events remain in the signal region while 8.4 ± 0.9 
events are expected from simulation. 
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3.2.4 Signal and background modelling 

For the dijet and lepton tag WW classes, the hypothesis testing is based on the number of ex- 
pected signal and backgrounds events only. Because of the impossibility to fully reconstruct 
the Higgs resonance and the small number of events expected in these classes, the limit extrac- 
tion is based on counting experiments in both the VBF and WH sub-channels. The number of 
expected signal events is evaluated from simulation, while the background contamination in 
the signal regions is estimated with methods based on data whenever possible. For the zero 
and one jet classes, the limit extraction is based on the shape of a multi-variate discriminant, 
optimised to maximise the difference between the signal and the WW background. The dis- 
criminant is built on the kinematics of the dilepton pair and the missing energy. For the signal 
case, the model is obtained from the simulation; for most of the backgrounds, the templates 
are taken from the simulation and cross-checked in control samples in data. For the W + jets 
background the nominal shape is derived from a data control sample. 

3.3 Diboson ZZ decay mode 

In the H — > ZZ^ — > 4£ channel [13] a search is made for a four-lepton mass peak over a small 
continuum background. The 4e, 4fi, 2e2}i sub-channels are analysed separately since there 
are differences in the four-lepton mass resolutions and the background rates arising from jets 
misidentified as leptons. The dominant irreducible background in this channel is from non- 
resonant ZZ production (with both Z bosons decaying to either 2e, or 2}i, or 2t with the taus 
decaying leptonically) and is estimated from simulation. The smaller reducible backgrounds 
with jets misidentified as leptons, e.g. Z + jets, are estimated from data. 

In the H — > ZZ — > 2£2v search [14j, events are selected by the presence of a lepton pair (ee 
or with invariant mass consistent with that of an on-shell Z boson, and a large E™ 1SS . 
A transverse invariant mass nij is defined from the dilepton momenta and E™ 1SS , assuming 
that E™ 1SS arises from a Z — > vv decay. A broad excess of events is searched for in the nij 
distribution. The non-resonant ZZ and WZ backgrounds are taken from simulation, while all 
other backgrounds are evaluated from control samples in data. 

In the H —¥ ZZ**) — > 2£2q search [15], events are selected with two leptons (ee or jiji) and two 
jets with zero, one, or two b-tags, thus defining a total of six exclusive final states. Requiring b- 
tagging improves the signal-to-background ratio. The two jets are required to form an invariant 
mass consistent with that of an on-shell Z boson. The aim is to search for a peak in the invariant 
mass distribution of the dilepton-dijet system, with the background rate and shape estimated 
using control regions in data. 

In the H — > ZZ — > 2£2t search [16|, one Z boson is required to be on-shell and to decay to 
a lepton pair (ee or jip). The other Z boson is required to decay through a tt pair to one of 
the four final-state signatures eji, er^, jir^, t^t^, where is a hadronically decaying t. Thus, 
eight exclusive sub-channels are defined. A broad excess is searched for in the distribution 
of the dilepton-ditau mass, constructed from the visible products of the tau decays, neglect- 
ing the effect of the accompanying neutrinos. The dominant background is non-resonant ZZ 
production whose rate is estimated from simulation. The main sub-leading backgrounds with 
jets misidentified as x leptons stem from Z + jets (including ZW) and top-quark events. These 
backgrounds are estimated from data. 

3.3.1 Signal and background modelling 

The limit calculation for the ZZ classes is based on the shape of the invariant mass distri- 
bution of the decay products, and the likelihood is written in terms of the estimated prob- 
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ability distribution function for the signal and the background. For H — > ZZ^ — > 4£ and 
H — > ZZS*> —> 2£2q, the signal shapes are described by means of analytical fits, based on the 
Crystal Ball function. For the H — > ZZ -4 2£2t class the likelihood is based on the binned distri- 
bution of the reconstructed visible mass. For the H — > ZZ — » 2£2v class the likelihood is based 
on the binned distribution of the transverse mass calculated using the visible decay products. 
The background shapes are extracted from data when possible, while for the irreducible ones, 
such as the electroweak ZZ or WZ production, the simulation is used. 

4 Results 

The statistical approach considered in evaluating the limit is the asymptotic CLg [32| with the 
profile likelihood ratio as a test-statistic [33). Given the narrowness of the Higgs mass peak in 
the 77 channel, which has a resolution approaching 1 GeV in the classes with the best resolu- 
tion, the search is carried out with steps of 0.5 GeV in the signal hypothesis mass in the range 
between 110 and 150 GeV. All known sources of systematic uncertainties are included in the 
likelihood model which is used for the limit setting. Systematic errors which are correlated 
between event classes (theory, luminosity, photon and trigger efficiency, etc) are included as 
common nuisance parameters. 

Following the prescription in [34|, the QCD scale uncertainties on the FP Higgs boson produc- 
tion cross section are increased with respect to those of the SM Higgs boson, to 5%, to cover the 
effects of electroweak corrections which have not yet been calculated. 

Figure [5] (left) shows the limit relative to the FP model expectation from the 77 sub-channels 
only, where the systematic uncertainties on the expected cross section and branching fraction 
are included in the limit setting procedure. The observed values are shown by a solid line. 
The contributions to the expected limit of each of the 77 sub-channels are also shown. The 
sensitivity of the search in this channel lies predominantly in the dijet tag sub-channel. The 77 
combined expected exclusion limit at 95% CL covers the mass range between 110-136.5 GeV, 
while the data exclude ranges from 110-124.5 GeV and 127-137.5 GeV. Figure [5] (right) shows 
the local p-value for 77 channel and each sub-channel, calculated from the asymptotic approx- 
imation [33 1, at 0.5 GeV intervals in the mass range 110-150 GeV. The local p-value quantifies 
the probability for the background to produce a fluctuation at least as large as observed, and 
assumes that the relative signal strength between the event classes follows the MC signal model 
for the FP Higgs boson. The local p-value corresponding to the largest upwards fluctuation of 
the observed limit, at 125 GeV, has been computed to be 3.6 x 10~ 3 {2.7a) in the asymptotic ap- 
proximation. When taking into consideration the look-elsewhere effect fl35| in the search range 
110— 150 GeV, the global significance of this deviation is 1.2 a. 

Figure|6]shows the 95% CL upper limit on the signal strength, a/ upp, of an FP Higgs boson, as a 
function of the Higgs boson mass, for the WW channel. The contributions from the individual 
sub-channels are indicated. The limit from the dijet-tagged sub-channel complements the 77 
search channels, excluding the FP hypothesis from 146 GeV to 196 GeV. 

The 32 sub-channels of the three decay modes, 77, WW, and ZZ, described in Section 3, are 
combined using the combination techniques described in Ref. l33l| to account for all statistical 
and systematic uncertainties and their correlations. The uncertainties consist of: theoretical un- 
certainties on the expected cross sections and acceptances for signal and background processes, 
experimental uncertainties in the modeling of the detector response (event reconstruction and 
selection efficiencies, energy scale and resolution), and statistical uncertainties associated with 
either ancillary measurements of backgrounds in control regions or selection efficiencies ob- 
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Figure 5: (Left) The 95% CL upper limit on the signal strength, a I o~pp, of an FP Higgs boson, 
as a function of the Higgs boson mass, for the 77 channel. The dashed line indicates the ex- 
pected median of results for the background-only hypothesis, while the two bands indicate the 
ranges that are expected to contain 68% and 95% of all observed excursions from the median, 
respectively. The asymptotic CLs method is used. Individual contributions to the expected 
limit for each of the channels are shown with dotted lines. (Right) The p-value of 77 channel. 
Contributions of individual sub-channels are also shown. 





Figure 6: The 95% CL upper limit on the signal strength, a/ c F p, of an FP Higgs boson, as a 
function of the Higgs boson mass, for the WW channel. (Left) Limit for the WW final state 
excluding the WH trilepton sub-channel. The contributions to the expected limit from the 
dijet-tagged and the single jet event classes are also shown. (Right) Limit for the WH trilepton 
sub-channel. The dashed line indicates the expected median of results for the background-only 
hypothesis, while the two bands indicate the ranges that are expected to contain 68% and 95% 
of all observed excursions from the median, respectively. The asymptotic CLg method is used. 
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tained using simulated events. 

The limit on the signal strength, a/ o~fp, of an FP Higgs boson, as a function of the Higgs boson 
mass, calculated using the asymptotic approximation, is shown in Fig. |7j together with the ex- 
pected and observed 95% CL limits for individual fermiophobic Higgs boson decay modes as 
well as for their combination. Checks at a few test points around 125 GeV have shown the calcu- 
lation to be consistent with values obtained by the full modified frequentist approach [36 1. The 
fermiophobic Higgs boson is excluded at 95% CL in the mass range 110-194 GeV. At 99% CL, 
we exclude the fermiophobic Higgs boson in the range 110-188 GeV, with the exception of two 
gaps: 124.5-127 GeV and 147.5-155 GeV. The sensitivity of the search lies predominantly in the 
77 channel below 140 GeV, and in the WW channel for the high mass search range. 

The local p-value as a function of the Higgs boson mass is obtained using the asymptotic ap- 
proximation for individual decay modes and for their combination, and is shown in Fig. [H] (left). 
For a few points, the p-value calculation is checked with the frequentist approach p3| and is 
shown to agree within the statistical error. The local p-value corresponding to the largest up- 
wards fluctuation of the observed limit, at 125 GeV, is computed to have a significance of 2.5 a. 
When taking into consideration the look-elsewhere effect [35J in the search range 110-300 GeV, 
the global significance of the deviation is 0.9 a. This deviation from the expected limit is too 
weak to be consistent with the fermiophobic Higgs boson signal, as can be seen in Fig.|8](right), 
which shows the observed signal strength for an FP Higgs, as obtained from the fit of signal 
plus background on data. In this fit the constraint on signal strength being non-negative is 
not applied, so that a negative value indicates an observation below the expectation from the 
background-only hypothesis. 




100 '150 200 250 300 100 150 200 250 300 

m H (GeV) m H (GeV) 



Figure 7: (Left) The observed and expected 95% CL upper limits on the signal strength, a I c F p, 
for the fermiophobic Higgs boson hypothesis as a function of the Higgs boson mass. (Right): 
The observed and expected 95% CL upper limits on the signal strength as a function of the 
fermiophobic Higgs boson mass for the three explored Higgs boson decay modes and their 
combination. 
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Figure 8: (Left) The observed and expected local p-values as a function of Higgs boson mass for 
the fermiophobic Higgs boson search. The local p-values for the three decay modes taken sep- 
arately, and their combination are obtained with the asymptotic approximation; the combined 
local p-value is validated at a few m\i (dots) using the full frequentist approach by generating 
ensembles of background-only pseudo-datasets. The dashed line indicates the expected com- 
bined local p-value, should a fermiophobic Higgs boson with a mass j«h exist. (Right) The 
best-fit value of the signal strength, a/ cr F p, for the combined sub-channels, as a function of the 
mass of a fermiophobic Higgs boson. 
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5 Summary 

Combined results are reported from searches for a fermiophobic Higgs boson in proton-proton 
collisions at -y/s = 7 TeV in its decay modes into vector bosons: 77, WW, and ZZ, where the 
Higgs boson production is restricted to vector boson fusion and associated production with a 
vector boson. The Higgs boson mass range explored is 110-300 GeV. The data analyzed corre- 
spond to an integrated luminosity of 4.9-5.1 fb _1 . The fermiophobic Higgs boson is excluded 
at 95% CL in the mass range 110-194 GeV and at 99% CL in the mass ranges 110-124.5 GeV, 
127-147.5 GeV, and 155-180 GeV. 
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